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Abstract 

The existence and the properties of isolated solutions to the relativis- 
tic Vlasov-Maxwell system with initial data on the backward hyperboloid 
t — — ^1 + |xP are investigated. Isolated solutions of Vlasov-Maxwell can 
be defined by the condition that the particle density is compactly sup- 
ported on the initial hyperboloid and by imposing the absence of incoming 
radiation on the electromagnetic field. Various consequences of the mass- 
energy conservation laws are derived by assuming the existence of smooth 
isolated solutions which match the inital data. In particular, it is shown 
that the mass-energy of isolated solutions on the backward hyperboloids 
and on the surfaces of constant proper time are preserved and equal, while 
the mass-energy on the forward hyperboloids is non-increasing and uni- 
formly bounded by the mass-energy on the initial hyperboloid. Moreover 
the global existence and uniqueness of classical solutions in the future of 
the initial surface is established for the one dimensional version of the 
system. 



1 Introduction 

Models representing isolated systems have an important role in physics. Even 
though no system can be rigorously considered isolated, it is common in physics 
to disregard in first approximation the effects that external agents, including 
the observer, have on the dynamics of the system under consideration. What 
exactly means that a physical system is isolated depends on the underlying phys- 
ical theory that is used to describe the dynamics of the model. In Newtonian 
gravity, for example, a compactly supported matter distribution is referred to as 
being isolated if the gravitational field energy at any fixed time decays to zero at 
large distances from the matter. This definition implies that the gravitational 
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potential generated by the mass distribution equals the convolution product of 
the mass density with the fundamental solution of the Poisson equation. How- 
ever if the dynamics of the same matter distribution is described using general 
relativity instead of Newtonian gravity, i.e., taking into account the relativistic 
effects, the precise definition of isolated sistem becomes more complicated and 
indeed is not yet fully understood, cf. |7|. 

In order to clarify the kind of difhculties that arise when we try to extend the 
concept of isolated system from pre-relativistic to relativistic physics, we study 
here the evolution of an isolated charge distribution in Maxwell's electrodynam- 
ics. As opposed to the electrostatic case, where the potential solves the Poisson 
equation and decays at the same rate in all directions, the finite, constant (in 
vacuum) speed of propagation of signals implies that a typical solution of the 
Maxwell equations will in general allow for a net flux of field energy coming 
from infinity and propagating along the null directions. The part of the field 
responsible for this energy flux is called incoming radiation. If the system is to 
be thought of as isolated, the incoming radiation should be ruled out by appro- 
priate boundary conditions. The precise form and relevance of these conditions 
will be discussed later. 

As regard to the matter model, we postulate a coUisionless plasma described 
by the relativistic Vlasov equation, which coupled to the Maxwell equations 
yields the relativistic Vlasov-Maxwell system. The reason for this choice is that 
solutions of the relativistic Vlasov-Maxwell system are in general more regular 
than for other matter models — such as perfect fluids — and this considerably 
simplifies the analysis in the sequel. In particular, the compact support property 
of the Vlasov matter field holds for all times provided it is required initially. 

1.1 The Vlasov-Maxwell system 

Consider a large ensemble of charged particles (plasma) which interact only 
by means of the electromagnetic forces that they generate collectively; exter- 
nal forces and collisions between the particles are assumed to be negligible. In 
the case of one species of particle, the average dynamics of the plasma is de- 
scribed through the particles distribution function f — f{t,x,p), which satisfies 
the Vlasov equation: 



Here (i,a;) G K. x M'^ denotes a Cartesian system of coordinates in the four di- 
mensional Minkowski space (i.e., a system of coordinates in which the Min- 
skowski metric has the canonical diagonal form) and p e R'^ is the momentum 
variable. The time coordinate t will be called proper time. {E,B) — {E,B){t,x) 
denotes the mean electromagnetic field generated by the plasma; the vector field 
E+px B is the Lorentz force and 



dtf+p-'7J+{E+pxB)-Vpf 
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denotes the relativistic velocity of a particle with momentum p. Units are 
chosen so that the mass and the charge of each particle, as well as the speed of 
light, have unit value. The electromagnetic field {E,B) is required to solve the 
Maxwell equations 

'dtE = S/xB-j, ^■E = p, 

dtB = -\/xE, V-B = 0, ^ ' 

where the charge density p and the current density j are given by 

p(t,x)^ / f{t,x,p)dp, j{t,x)^ pf{t,x,p)dp. (1.4) 



The set of equations Hl.l|l - H1.4|) is the (relativistic) Vlasov-Maxwell system. 
Sufficiently regular solutions of this sistem satisfy the continuity equation 



as well as the energy identity 



where 



dtp + V-] = 0, (1.5) 
dte + V-p = 0, (1.6) 



e{t,x) = J ^^TWfdp+^\E\^ + ^\Bf^ p{t,x) = Jpfdp + ExB, 

are the local energy and the local momentum, respectively. Due to its relativis- 
tic character, the Vlasov-Maxwell system allows for the presence of radiation 
fields. The radiation is defined as the part of the electromagnetic field which 
carries energy to null infinity. It is distinguished in outgoing radiation, which 
propagates energy to future null infinity, and incoming radiation, which propa- 
gates energy to past null infinity. As for any solution {f{t,x,p),E(t,x),B{t,x)) 
of Vlasov-Maxwell, the triple {f{~t,x,—p),E{—t,x),—B{—t,x)) gives a new so- 
lution thereof (time reflection symmetry), solutions of this system contain in 
general outgoing as well as incoming radiation. Our main interest is in the so- 
lutions of the Vlasov-Maxwell system which describe an isolated plasma. The 
precise definition of this class of solutions will be given in Section |21 The phys- 
ical idea behind the concept of isolated solution is that the system is not hit 
by energy coming from infinity. This condition can be achieved by requiring 
that (i) the particle density is supported away from past null infinity and (ii) 
there is no contribution from the incoming radiation in the electromagnetic field 
{no-incoming radiation condition). Note that since radiation propagates energy- 
momentum at the speed of light, the correct set of infinity where to impose the 
condition for having an isolated system is past null infinity. 



1.2 The initial value problem with data on a backward 
hyperboloid 

As the definition of isolated systems is related to the behavior of the solution 
as i— >— oo, the standard framework of the Cauchy problem with data at t = 
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does not seem to provide a suitable mathematical setting for the analysis of 
such solutions (local isolated solutions of the Cauchy problem are meaningless) . 
In the framework of the initial value problem, it is more natural to assign the 
initial data of isolated solutions on a surface which cuts past null infinity. In 
this surface was chosen to be the past light cone with vertex in the origin 
of the Cartesian system of coordinates. This choice leads naturally to consider 
the foliation of Minkowski space by the family of past light cones t—\x\=v 
and to study the Vlasov-Maxwell system in the advanced time coordinate v. 
However this situation is complicated by the fact that, since the coefficients 
of the Vlasov equation are discontinuous at the vertex of the cones, smooth 
solutions of the Vlasov-Maxwell system in the advanced time coordinate do 
not exist in general. Considering weak solutions to the problem does not seem 
completely satisfactory, because the singularity at the vertex of the cones is only 
an artefact due to a bad choice of coordinates and not a "real" singularity. 

Motivated by the previous remarks, we continue the analysis on isolated 
solutions of Vlasov-Maxwell started in [H] by employing a different approach, 
namely by studying the initial value problem with data on the backward hyper- 
boloid t — — jxp. Besides solving the problem on the regularity of isolated 
solutions mentioned above, the present paper also contains several results which 
were left open in jHj or merely conjectured therein. 

Let (/n,i?ni-Sn) denote the restriction of the solutions of Vlasov-Maxwell 
on the backward hyperboloids t = T— ^1 + \x\'^, r^O. In order to make trans- 
parent the relation between (/n,i?n:-Sn) and {f,E,B), the coordinates on the 
hyperboloids will be chosen by orthogonal projection on the hypersurfaces of 
constant proper time. Hence 



fn {T,x,p)^f{T- ^/T+W, ^,P), {En,Bn){r,x)^{E,B){T- ^/T+W, ^) ■ 



In terms of these new dynamical variables, the Vlasov-Maxwell system takes 



the form 



{l+p-x)drfn+P-'^xfn + {En+P^Bn)-S/pfn^O, 



(1.7) 



dr{En -XX Bn) = V x Bp - jn, 
dr{Bn+xxEn) = -^xEn, 



(1.8) 
(1.9) 



dr{En-x)+\/-En^Pn, 
driBn-x) + V-Bn^O, 



(1.10) 
(1.11) 



where 




(1.12) 



Initial data are given at r = and denoted by 



f^^ix,p)^fn{0,x,p), El^\x)^En{0,x), B'^{x) = Bn{0,x). 
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Later we shall discuss the equivalence of the system above with the evolution 
equations (|1.7() - (|1.9() and a set of constraint equations on the initial data. We are 
interested in the question of existence and uniqueness of classical solutions in the 
future, i.e., for rG [0,oo), which match the initial data at r = 0. These solutions 
will be referred to as "global" , although they are defined only for r > (note that 
p.7|l - H1.12|l is not time symmetric) . Obviously, data on a backward hyperboloid 
are in general not sufficient to determine a unique solution in the future, since 
the intersection between the initial surface and the domain of dependence of the 
solutions on a space-time point sufficiently far in the future is not a compact set 
(precisely, uniqueness is lost in the region of Minkowski space in the future of 
the past light cone t = — |a;|). Based on causality arguments, it is expected that 
a unique solution of H1.7|l - H1.12|l is obtained by supplying the initial conditions 
with data at past null infinity (i.e., at |a;|^+oo on the surfaces t — const.). 
That this is the case will be proved in Section |2| for the one dimensional version 
of the system. In the class of solutions which admit a continuous extension at 
past null infinity, the no-incoming radiation condition is equivalent to giving 
zero data at past null infinity. 

1.3 Outline of the paper and outlook 

The outline of the present article is as follows. In Section|2]various consequences 
of the local conservation laws for H1.7|l - H1.12|l and of the no-incoming radiation 
condition are collected. The results of Section |21 are conditional, as they assume 
the existence of global classical solutions, and improve those of [B] in several 
aspects. In particular, the boundedness of the mass-energy on the forward 
hyperboloids, which was left open in is established here by integrating the 
identities H1.5(l - H1.6() in the interior of suitable space-time regions. Moreover we 
prove that the energy of isolated solutions on the surfaces of constant proper 
time is preserved and equals the energy on the initial hyperboloid, a property 
which was only conjectured in [S]. In Section f2.6l we state a global existence and 
uniqueness theorem for spherically symmetric solutions. The proof is based on 
classical arguments introduced in ,1. ,11 and is therefore omitted. Note that in 
spherical symmetry the magnetic field vanishes identically (if decay at infinity 
is imposed) and the Maxwell equations reduce to the Poisson equation for the 
electric field. Hence there is neither incoming nor outgoing radiation in spherical 
symmetry. In order to introduce radiation effects, we consider in Section 13 the 
Vlasov-Maxwell system in one spatial dimension and two dimensions in velocity 
(the so-called "one and one-half dimensional system). In this case we are able 
to prove the existence and uniqueness of a global classical solution with given 
data at r = and data at past null infinity. Moreover we show that if and only 
if the data at past null infinity vanish, the solution satisfies the no-incoming 
radiation condition. 

To conclude the Introduction, we review some important results on the 
Vlasov-Maxwell system. Existence and uniqueness of local solutions to the 
Cauchy problem has been proved in In it is shown that a classical 

solution exists in any interval of time in which the momentum support of the 
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particles density is bounded. Different proofs of this result are given and [HI- 
In JH| it is proved that such a continuation principle of local solutions hold for 
a wide variety of linear equations coupled (in a non-linear way) to the Vlasov 
equation. The result in was applied to prove global existence and unique- 
ness under different smallness assumptions on the initial data (cf. |12[ 1151 [TO] ) 
and for arbitrarily large data in one and two space dimensions [131 114| . Exis- 
tence, but not uniqueness, of global weak solutions in three spatial dimensions 
is proved in 

Isolated solutions of the Vlasov-Maxwell system were first studied in ^ 
using the retarded solution of the equations. For related results based on 
approximation methods see |2] . 

2 General properties of classical solutions 

In this section we study several properties of global solutions of H1.7|l - H1.12|l 
satisfying the regularity condition 

fneC\[0,oo)xR^xR^), En.Bn eC\[0,oo) xR^) (2.1) 

and so they are solutions of H1.7|) - (|1.12|l in a classical sense. It is also assumed 
that /n has bounded support in the momentum, precisely 

Vnir) = sup{ v/1 + bp : fn{s,x,p) 7^ 0, s r, x G R^} < cx), VrG [0,oo). 

(2.2) 

In particular, all the integrals in the momentum variable in the sequel are un- 
derstood to be extended over a (time dependent) compact set. For the sake of 
reference we define 

Coo = {{fn,En,Bn) which satisfy (EHJ-lESl}- 

The precise mathematical definition of isolated solution of the Vlasov-Maxwell 
system can now be given. We denote by S'r = {2; G M"^ : = r} the sphere with 
center in the origin and radius r>0, by dSr the invariant volume element 
thereon and by k = x/\x\ the outward unit normal to the sphere. 

Definition 1 A solution (/n,i?n,-Bn) gCoo of \l.'T\j - ^.l^) is said to be isolated 
if 

(A) There exists a constant i?o > such that 



pn(r,x) = 0, for ^1 + \x\^;^Ro + -t; 

(B) The no-incoming radiation condition (NIRC) in the future is satisfied, 
namely 
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Remark 1 A local isolated solution in the interval [0,T), T> 0, can be defined 
by imposing the conditions (A) and (B) only for t,ti,T2 £ [0,T). 

Remark 2 The limit in the left hand side of (|2.3|) defines the energy carried 
out by the electromagnetic field to past null infinity in the interval [Ti,r2] of 
the time coordinate r. As we are interested only in going forward in time, we 
require the NIRC only in the future of the initial hyperboloid t — 0. 

Remark 3 The condition (A) means geometrically that the matter is supported 
away from the region of past null infinity in the future of the initial hyperboloid 
(however the matter support may intersect future null infinity). Although this 
assumption could be relaxed and replaced by a decay condition that implies the 
absence of kinetic energy at past null infinity, we shall avoid to do so, as we feel 
that such a condition does not add any phyisically interesting new feature to 
the concept of isolated system. 

It turns out that the condition (A) in Definition is satisfied by any solution in 
the class Coo if the initial datum /p" is chosen with compact support in x and 

defining 

Ro = inf {r > : f^\x,p) = 0, ^/T+W ^ r, p £ R^}, (2.4) 

see LemmaEl For this reason, the initial data are throughout assumed to satisfy 

Oi^f^'eCliR^xR^), E'^,B'^'eC^{R^). (2.5) 

At this point it is worth to emphasize that even if the particles density were 
given with compact support on the surface i = (as it is usually the case when 
studying the Cauchy problem), there is no guarantee a priori that it will have 
compact support also on the backward hyperboloid t = — -^/l + \x\'^. This already 
indicates that the class of solutions of H1.7|l - H1.12|l with data as in H2.5() is strictly 
contained in the class of solutions with data prescribed at t — 0. 

Remark 4 Definition can be adopted also for other matter models coupled 
to the Maxwell system, such as perfect fluids described by the Euler equations. 
However in the latter case it is not enough to give a charge density with compact 
support in the initial hyperboloid to make the condition (A) satisfied. To this 
purpose one needs to formulate an initial- (free) boundary value problem for the 
Maxwell-Euler equations, for which there is no existence result available. 

We split our analysis in six different subsections. The results of Sections l2. H and 
12. 61 make use explicitely of the Vlasov equation, whereas the content of Sections 
I2.2H2.5I applies to any matter model for which the local conservation laws of 
mass and energy hold. 

2.1 The Vlasov equation 

Let us start by pointing out some basic properties of /p. The estimate 
-1,^^ 1 , P-x \p\\x\ 

1+p-X = 1+ , , ,^ , , ,o 
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bl 



shows that for solutions in the class Coo i the equation (|1.7|l is equivalent to 



drfn + — ■ V,/n + — Vl + \p^En+px Bn ■ Vpfn = 0, (2.6) 
Pq Po ^ J 



where po is defined by 



The characteristics of the differential operator in the left hand side of (|2.6|l 
are the solutions of 



p^ — (^/l + \pVEn+pxBr 



1 



Po 



Po 



(2.7) 



and we denote by (X,P)(s) the characteristic curve satisfying {X,P){t) = (x,p); 
note that {X,P) also depends on {t,x,p), but this is not reflected in our notation. 
Let F={Fi,F2) denote the right hand side of (|2.7|) . that is 



Fiix,p)^^, F2{x,p) = — Ul + |p|2Sn +P X B, 
Po Po ^ 



(2.8) 



Since F is a function when {Ed,Bd) G C^, then the characteristics are smooth 
and one can write the solution of 12.6|l as 



/n(T,x,p) = /^"((X,P)(0)). 



(2.9) 



In particular, /n remains non-negative for all times and ||/n(''')||oo ^ ||/n"||oo- 

We shall now derive the conservation laws satisfled by the solutions of H2.6|l . 
A straightforward computation shows that the vector (|2.8|l satisfies 



(2.10) 



[V(,,p) • F] {s,X{s),P{s)) - -^log {l + P{s) -Xis) 



where P — + |Pp. In fact, each side of H2.1()|l equals, along characteristics. 



{i+p-xy 



\p\^-{p-x)^ 



From H2.10|l we deduce 

-d{X,P){s) 



/l + \ 12 ■ " ^^P^ - (P X ^) ■ -Bn) 



det 



d{x,p) 



1 +p- X 



1 + P{s)-X{s) 



Hence using H2.9|l the next lemma follows. 
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Lemma 1 For any solution of in the class Coo and for any 

surable function Q:R— s-M, we have 



(3(/n)(l +p- x) dpdx = const. 

In particular, by choosing Q{z) = z'^ , q ^ 1, 

/n(T)||L, (R3xR3) =const. (2.11) 
In the next lemma we estimate the a;— support of fa- 

Lemma 2 Let (/n , i ) G Coo be a solution of ji.7l )- ji. with initial data 
fr3|] . Then 

fn{T,x,p)^0, for ^l + \x\^^Ro + ^T, r^O, 
where Rq>0 is given by i2.4\ )- 

Proof: For all ^ s ^ r we have, by the first equation in H2.7|l . 

' Js ^1 + \P{t')\^+P{t')-X{t') 

ii+ ^1 + \P{t')\^+P{t')-X{t') 



where X = X/^1 + |X|2 andX+ = {r'e [0,r] : (P-X)(t') ^0}. Since y/T+W> 
\p-x\ we obtain 



^ ^1 + \X{S)\^ + -{T-S). 



For s = this implies \J1 + |a;p ^Ri:i + \t in the support of /n . □ 

Corollary 1 A solution (/n , i?n j ) £ Coo with initial data \2.5]) is isolated if 
and only if it satisfies the NIRC in the future. 

2.2 The Maxwell equations 

Next the Maxwell equations (|1.8I) - (|1.11|) will be considered. A basic fact is that 
they are equivalent to a set of constraint equations on the initial data and a set 
of evolution equations. As opposed to the case of initial data on a past light 
cone considered in 0, data on a backward hyperboloid are restricted by the 
same number of constraints as for the case of data given on the surface t = 
(in [S] there appear additional constraints because the initial data are given on 
a characteristic surface). We consider the system (|1.8|l - (|l.ll|) for a given pair 
(pn,in) G satisfying the equation 

5r(pn+jn-S) = -V-in. (2.12) 
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Scalar multiplying p.8|l - p.9|l by x and using p.lO|) - (|l.ll|) we obtain 

V ■En = ipn+jn-x)-x-V X Bn, V ■ B^^x -V x E^. (2.13) 

A simple vector algebra computation shows that the equations (|2.13|) are sat- 
isfied for all times provided they hold at t = and the functions i?n,5n;Pn,jn 
satisfy (|1.8|) . H1.9|) and H2.12|l . Hence we have the following 

Proposition 1 Let {pr],jn)^C^ be given which satisfy 1^2. IS)) . Then {En,Bn) 
is a solution of with initial data {E™,B^f^) if and only if it is a 

solution of H1.8p - I^l.y\} with initial data satisfying the constraint equations 

V-El^^ip^+f^^-x)~x-VxBl^, V-Bi^^ = x-VxEi^\ (2.14) 

where (Pn > Jn") ~ (Pn,in)|r=o- Moreover the solution satisfies \2.1!^) for all t^O. 

Proposition n] applies to solutions of H1.7|l - H1.12|l with p"' ~ j fn'^P j"^ = 
jpf™dp, since the validity of (|2.12|) follows by the Vlasov equation (or by a 
simple change of variable in p.5|l ). 

2.3 Evolution of the mass functions 

Our next purpose is to study the evolution of certain mass functions associated 
to a solution in the class Coo- For t^O fixed and (5 = 0,1,2, consider the hy- 
persurfaces in Minkowski space which, in the Cartesian coordinates (i,a;), are 
defined by 

S5 (r) = { (i , x) : i + (1 - (5) VT+N^ = r, X e ^ 

It follows that S5(t) is a backward hyperboloid for (5 = 0, a surface of constant 
proper time for 5=1 and a forward hyperboloid for 5 = 2. The notation [/i]s_5(r) 
will be used to indicate that the function h = h{t,x) has to be evaluated on the 
surface T,s{t). For r ^ and r > fixed we set 

nsiT,r)= / [p+{l-5)j-x]^^(^r)dx, 

J Br 

where B^ denotes the ball with center in the origin and radius r > 0. In terms 
of a solution (/n,£^n,^n) gCqo of H1.7|) - (|1.12|l . the functions ug are given by 

ns{T,r) = J [pn + il-5)jn-x] (t + 5 ^/l+Jx^ , x^ dx. 

We set 

Ns{t) = lim ns{T,r), 

r — >oo 

which is the total mass of particles on the surface ^^(t). The limit in the 
previous definition exists, since ^^(r, •) is non-decreasing. It is clear that A'o can 
be defined also for local solutions, whereas A^i and N2 make sense only for global 
solutions. Let A^q" = A^o(0), which depends only on the given initial datum /™ 
and is bounded for data as in (|2.5|l . 
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Proposition 2 Let {fn,Er,,Bn) GCoc be a solution of j i. 7| )- i) i. with initial 
data 1^2. 5\) . Then the following holds: 

(1) iVi (r) = iVo(r) = N'^, for all r ^ 0; 

N2{t2) ^ A^2(ti) sC TV'", for all «C n < Ta; 

l',?; //sup{7'n(T),re [0,c»)}<oo, then N2{t) = N'^ , for allr^O. 

Proof: The equality No{t) = iV™ corresponds to (|2.11|) for q — 1. Integrating the 
continuity equation (|1.5|) in the interior of the space-time region 



7^l (r, r) = { (t , x) : T - vT+N^ s; t T, X e 

and applying the divergence theorem we obtain 



ni(T,r) =no(T,r)^ 
By Lemma El and fTT^ . 



T+Vl+r 



[in ■k]{T' ,x)dSrdT' . 



(2.15) 



ni{T,r) —no{T,r), for r > (2i?o + T) — 1 



1/2 



which in the limit r—Kx implies A^i (r) = A^o (t) . By H1.5() , the function n2 
satisfies the equation 



dr 



-dr 712 = - 



[ph2iT)dSr 



(2.16) 



Integrating along the characteristics of the operator in the left hand side of 
H2.16|l we obtain 



"2 T2 



iVl + r^-T2f-l 



1/2 



1/2 



for T2 ^ Ti and r ^ ■\/ (1 +T2)^ — 1. Letting r — ^ oo in the previous inequality 
shows that N2{t) is non- increasing. To complete the proof of (2) it remains to 
show that A^2(0) ^N™. To this purpose we integrate (|1.5|l in the interior of 



7^2(r) ^{{t,x):Qi^t^ VT+W, ^ ^ Br-t}, r>l, 
and apply again the divergence theorem to obtain the identity 

r^-l' 



712 



:ni(0,r)- 



2?- / .10 Js,^-, 



[p + j-k] {t,x)dSr-tdt. 



Since the integral in the right hand side is positive, we obtain the inequality 
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2r 



<ni(0,r), 
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whence, in the limit r^oo, N2{0) ^ Ni{0) = . To prove (3) we use the 
identity 



n2{T,r) = ni{T,r)- 



+ VT+7' 



[jn-k] {t',x) dSrdr', 



(2.17) 



which foUows integrating H1.5|l in the interior of the space-time region 



7^3 (r, r) = { (t, x) : T r + yi+N^, X e B J . 

Let D = sup{'PniT),Te[0,oo)} and x eR^ : fn{T,x,p) ^0 for somcpeM^. The 
assumption D <oo imphes 



V 1 + \p\ ^ — — \p\ ^ — — (p-x) 

in the support of /n and so, as in the proof of Lemma |21 

P(r')-X(r') 



1+ ^1 + \P{t')\^+P{t')-X{t') 

T, 



-dr' 



for aU (a;,p) esupp/n(T), where X+ = {T'e [0,r] : {P ■ X){t') > 0}. This imphes 
that 



fn{T,x,p)^0, for ^l + |a;|2^iT!o + ir, a=— S=== 

^ D + \/\ + 



°4 



Whence 



[jn-fc] (T',a;) dS^dr'^O for r> 



Ro + ar 
l-2a 



-1 



1/2 



□ 



Thus (3) follows letting r~*(X) in (|TT7|) . 

If N2{t2) < N2{ti), for ^ Ti < r2, the difference iV2(Ti) - ^^2(7-2) defines the 
mass (or number) of particles lost at future null infinity in the interval [Ti,r2]. 
This outgoing mass is measured by observers at future null infinity. Under the 
assumption in (3) of Proposition |21 no particles reach future null infinity and 
therefore the outgoing mass is zero. We remark that a uniform bound on VniT) 
seems to hold only under severe assumptions, such us spherical symmetry (see 
Section ing , or small data, see p iTHllTnj . 

Remark 5 The conclusions of Proposition |21 are valid even if the solution does 
not satisfy the NIRC. 
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2.4 Evolution of the energy functions 

Now we define the energy of a solution of the Vlasov-Maxwell system on the 
hypersurface ^^(t) as 

Ms{t) ^ lim ms{T,r), ms{T,r)^ [e+{l-S)p-x]^^,r)dx, 

and we set Mq^ — Mq{0). The next goal is to study the evolution of these 
energy functions. We proceed along the same lines as for the evolution of the 
mass functions. 

Proposition 3 Let (/n,£'n,-Sn) gCoo be an isolated solution of il/I}) - !^!.!^} 

(or, equivalentely, let the solution satisfy the NIRC in the future and the initial 
conditions ) and assume M^<oo. Then the following holds: 

(1 ) Ml (t) = Afo(T) = M'^, for all r ^ 0; 

(2) M2{t2) < M2(ti) ^ A4^, for all < n T2. 
Proof: By we have 

drmo = - [pn-k]{T,x)dSr- 

JSr- 

Integrating in the interval [0,r] we obtain 

mo(T,r) = mo(0,r)- / / [pn • fc] (T',a;)d5,dT'. (2.18) 

"'0 JSr 

In the limit r—^oo, the integral in the right hand side converges to zero by the 
support property of /n and the NIRC. This shows that Mo(r) is constant. Now, 
again by H1.6|) . 



{dr±dr)mi^ J [±e~p-k]^^^^~^dSr 



Since the right hand side of the previous equation is positive in the plus case 
and negative in the minus case, it follows (integrating along characteristics and 
letting r^oo) that Mi is constant. To complete the proof of (1), it remains to 
shows that Mi(0) = M™. Integrating the energy identity H1.6|) in the interior of 
the space-time region 



7^4(r) = {{t,x) : -y/l + W^t s: 0, X e Br+t}, r>l, 
and applying the divergence theorem we obtain 

mi(0,r) = mo fo,^— + / / [e-p ■ k]{t,x)dSr+tdt 



2r 



+ t 



^mo(0,^— ). (2.19) 
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On the other hand, integrating (|l.t)|) in the interior of the space-time region 
7^5(r-i,r2) = {(i,x):-V'TTN^=^<<min(0,ri-A/TTN^),xe-B^2}, 



ri>l,r2> \Jrf-l, 

we obtain 

mi ^0,y'rf~T^ = mo(0,r2)-^ J [pn • fc] (r, x) d^^^ ^^r 

[en+pn-x]{ri,x)dx 

As we are considering isolated solutions, the integral in the right hand side 
of the latter inequality converges to zero in the limit r2— s-oo, which entails 
TOi(0,r) ^Mq", for aU r >0. Combining with H2.19|l we get 

mo(^0,^^^ ^mi(0,r)<AC, Vr>l, 

whence Mi(0) = M^". The equation 



-r 



\^dr^^^drjm2 = -^^ NE2(r)rf5r, (2.20) 

shows that A/2 is non-increasing. The estimate Al2{0) ^ A/j", which concludes 
the proof of (2), follows by the same argument used in Proposition[21to establish 
the bound A^2(0) ^ A'^", i.e., integrating the energy identity (|l.t)|) in the interior 
of the space-time region 7?.2(t). □ 

If M2(t2) < Ar2('''i), for ^ Ti < T2, the difference Ai^2('''i) — Af2('''2) measures 
the energy dissipated at future null infinity in the interval [Ti,r2]. This outgoing 
energy is measured by observers at future null infinity. 

Remark 6 One clearly expects that the absence of outgoing energy implies the 
absence of outgoing mass. To see this, note that H2.16|l and H2.20|l imply 



\^T- " ' drj (m2~n2) = - " " y J ([e]E2(r)-[p]s2(r)) dSrS^O, 

which entails Af2(T2)- Afi(Ti) s^iV2(T2)- A^i(ti) <0. Hence A/i(ti) Af2(T2) 
implies A^i (n ) = A^2 (t2 ) . 
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Remark 7 If the solution is not isolated, then the energy on the backward 
hypcrboloids is non-decreasing. For the function mg satisfies the equation 



dr-i drjmo = J [e\^g(^r)dSr^O, 

whence Mo{t2) ^ Mo{ti), for all T2^ti. 



2.5 A uniqueness theorem 

Let us show that the NIRC guarantees uniqueness of solutions to the Maxwell 
equations with data on a backward hyperboloid. 

Proposition 4 Let pn^jn G be given such that the continuity equation 1^2.1'^) 
is satisfied. Let G C^(R'^) be given such that the constraint equations 

\2.1^^ are satisfied and Mo(0) <<X), where 

Mo(r) = i j (|SnP + |SnP + 2(SnXBn)-^)rfa:. 

There exists at most one solution {Er\,Br\) of the Maxwell system 

E221) satisfying D and {Er^,Bn)\t=o^{.E]^,B\^)- 

Proof: The difference {6Ed,SBd) of two solutions with the same data satisfies 
p.8|l - p.ll|l with Pd =jn = and zero data; by Proposition^ {6Ed,SBd) solves 
the equations 

dr{SEn-xxSBn) = '^xSBn, dr{SBn+x x dEn) = -S/ x 5En, (2.21) 

V-6En^~x-Vx6Bn, V ■ 6Bn=x-V x SE^. (2.22) 

On the other hand, the energy of {6Ed,SBd) on the backward hypcrboloids is 
zero; this implies 

^ \SEn-x\'^ + \SBr,-x\'^ + \SEr,-xxSBn\^ + \5Br,+xxSEn\'^ 

= {\6En\' + \SBn\')+A{SEnxSBn)-x 

l + \x\'^ 

i^i(^^\SEn\^ + ^\SBn\^ + {SEnx6Bn)-x^=0, 

whence dEr]~xx SBr]=SBr)+xxSEr]—0, which together with H2.21|l - H2.22|l 
yields SEn^SBn = 0. □ 



2.6 Spherically symmetric solutions 

In spherical symmetry we have V x Ed = V x Bd — and so, by the second 
equation in 12.13|l . V-i?n=0. Under the additional boundary condition 
lim|2,|^oo^n =0, this implies that the magnetic field vanishes identically. Thus 
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all spherically symmetric solutions satisfy the NIRC and so, choosing f™ of com- 
pact support, they are all isolated. Moreover, by the first equation in H2.13|l . 

k r 

En{T,x) = — / ipn+jn-x)iTyy^dr' 

{pn+jn-x){T,y)dy, (2.1) 



r- 

1 f (x-y) 



Att J \x — y 



the second equality being valid in spherical symmetry. By abuse of notation 
we use the same symbol to denote a spherically symmetric function in spherical 
and Cartesian coordinates. The Vlasov equation reduces to 

gr /n + — ■ V^/n + ^ En ■ Vpfn = 0. (2.2) 
Po Pq 

Note also the conservation of angular momentum; along characteristics, 

^\xxp\^^0. 
as 

In the spherically symmetric case we have the following global existence theorem. 



Theorem 1 Let f™ eC^{M.^ xM.'^) be spherically symmetric; there exists a 
unique (spherically symmetric) fn & C-^([0,oo) x R'^ x MP) solution of f2.1}) - f2.^ 
such that fr]{0,x,p) = f™(x,p). Moreover, there exists a constant D>0, depend- 
ing only on bounds on the initial datum, such that 

VniT)^D. (2.3) 

The proof of Theorem ^ is a straightforward adaptation of the argument given 
in [111 Theorem 11] . Since the details of this adaptation have already been given 
in [H], they will not be repeated here. Note however that in |H] an extra condi- 
tion on the initial datum is assumed, namely that there are no particles with 
arbitrarily small angular momentum at time zero. By the conservation of an- 
gular momentum, this condition is preserved in time, and since the maximal 
momentum of the particles is uniformly bounded, the conclusion follows that 
the particle density vanishes in a neighbourhood of the axis r = 0. This property 
was necessary in the analysis of ^6^ because, as we mentioned in the Introduc- 
tion, the particle density is singular on the vertex of the backward light cones. 
In the case considered here, however, this problem does not arise, and no further 
condition is needed on the initial datum. 

Corollary 2 For the solution of Theorem^ 

N2{T)^N,iT) = NoiT) = Ni;\ 
M2 (t) = Ml (t) = Mo (r) = A^" . 
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Proof: The equalities iVi (r) = TVq (t) = iV^" and Mi{t) ^ Mo{t) = M'^" hold for 
all isolated solutions, whereas N2{t) = iVg" follows by (|2.3|l . see (3) of Proposition 
13 For B = Q, the function m2 satisfies an equation similar to H2.17f) . namely 

"^2(T,r) = ml(r,r)- / / p-kfnir' ,x,p)dpdSrdT' . 

As in the proof of (3) of Proposition |21 the uniform estimate (|2.3|l implies that 
the integral in the right hand side vanishes for r large enough, whence the 
identity M2 (r) = Mi (r) follows by letting r — > cx) in the previous equation. □ 



Remark 8 From a physics point of view, it would be desirable that the emis- 
sion of outgoing energy is a generic property of isolated solutions. We conjec- 
ture that spherically symmetric solutions are the only isolated solutions of the 
Vlasov-Maxwell system for which the energy M2(r) is constant. 

3 The one dimensional system 

In this section we consider the Vlasov-Maxwell system in one spatial dimension 
and two dimensions in velocity, which is also known as the "one and one-half 
dimensional Vlasov-Maxwell system" . We remark that, besides greatly reducing 
the mathematical complexity of the problem, one dimensional Vlasov-Maxwell 
systems have applications in plasma physics — e.g., to model laser-plasma in- 
teractions or electron beams — and in numerical simulations, see jl(J| and the 
references therein. In Cartesian coordinates, the Vlasov-Maxwell system in one 
spatial dimension and two dimensions in velocity takes the form 

dtf+pAf+ [El + Bp2)dp, f+iE2- Bpi)dp, / = 0, (3.1) 

dtEi^~ji, d^Ei^p-n, (3.2) 

dtE2+d^B = -j2, dtB + d,E2 = 0, (3.3) 

ipJij2)^ J{l,Pi,P2)fdp, (3.4) 

see |13| . The particle density in phase-space is the function f ~ f{t,x,p), where 
{t,x)GM.'^ denotes a Cartesian system of coordinates in the two-dimensional 
Minkowski space, p—{pi,p2)&^'^ is the particles momentum and, for all zSM, 
we set 

In particular, p= (^1,^2) is the relativistic velocity. The electric field is given by 
E= (£'i,i?2) = E{t,x), while the magnetic field is the scalar function B — B{t,x). 
The function n = n(x) is a background density that will be chosen in such a 
way that the component Ei of the electric field vanishes at infinity (the latter 
condition is necessary for the existence of solutions with finite energy). The 
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continuity equation and the energy identity for the system (|3.1|I - H3.4|I take the 
form 

dtP + dxji^O, dte + d^p^O 

where 

eit,x) = ^\E\^ + ^B^ + J V^TWfdp, p{t,x)=E2B + Jpifdp. 

As in the three dimensional case, we start by reformulating the system 
H3.3|l in terms of new dynamical variables evaluated on the backward hyper- 
boloids t + vl + a? = r, so that the limit at past null infinity becomes x^±oo; 
define 



fn{T,x,p) = f{T-Vl + x^,x,p), {pn,jn){T,x)^{p,j){T- Vl + x^,x), 

Un{T,x)^Ei{T^VT+^,x), (0n,V'n) = ^ {E2 + B,E2-B) [t- ^/l + x^,x). 
The system H3.1|l - (|3.3|l in terms of these new unknowns reads 

{l+pix)drfn +Pidxfn + [Un + {,(t>n ~ ikn)p2]dpjn 

+ [0n(l-pi)+V'n(l+Pi)]5p,/n=O, (3.5) 

drUn^-Jn, dxUn=Pn+j^x-n, (3.6) 

dr + T^dx)4>=--^ (dr--Kdx]ij=--^ (3.7) 

l + x J 2[l + x) \ l — x ) 2(1— .t) 

The NIRC in the future for the system (|X^ - l|X7jl is 

Ymi l\[,i>l-^l)(T,-T)~{<i>l-^l){T,r)\dx^Q Vri,r2>0. (3.8) 



The main result of this section is a global (in the future) existence and unique- 
ness theorem of classical solutions to the system H3.5ll - H3.7(l . Moreover a neces- 
sary and sufficient condition on the data of the problem is given such that the 
solutions satisfy (|3.8I) . In order to give a precise formulation of our results, we 
need first to introduce some preliminaries and notation. Firstly we shall assume 
throughout that the background density n{x) is smooth with compact support 
and is neutralizing, namely 

[(p + jnx)(0,a;)-n(x)] dx = 0. (3.9) 

This is the only possibility which leads to a finite energy solution, see Remark 
El. Let (en, pn )(''', a:) = (e,p)(T — Vl + a;^,^); the energy identity takes the form 

5r(en+pn^) = -5:.pn. (3.10) 
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The mass-energy on the backward hyperboloids in one spatial dimension are 
given by 

Nq{t)= Mm no{T,r), no{T,r)^ [pn+jnx){T,x)dx, 



Mo(r) = hm Too (T,r), mo{T,r)= {en+pnx)iT,x)dx 
and as before we set iVo(O) = N^"" and A'/o(0) = A/^". Note that 



= / lpofndxdp + J 



iL/^(0,x) + (C)'(l + S) + (C)'(l-^) 



dx, 



where we denoted 

Po = V^TW+Pi^, = (/n,0n,^n)|..o- 

For T>0, let Ct denote the set of quadruples (/n,[/n,0n,^n) of functions 

/n:[0,r)xMxR2->[0,cx)), Un,(j>n,iJn-[0,T)xR~^R, 
with regularity /n , C^n , 0n , "V^n S and such that 



Vn (r) = sup{ v/l + b|2 : f{s,x,p) ^ 0, s G [0,r] , x G K} < oo, Vre[0,r). (3.11) 

The solutions of (|3.5|l - (|3.7|l that we construct belong to the class Coo- This set 
of solutions is the analogue of what we used in the three dimensional case. In 
particular, for solutions in the class Coo the results of Section El hold with the 
obvious modifications in the one dimensional case. In order to obtain a unique 
solution in Coo, an initial data set has to be prescribed. We assume 

(I) : M X ]R2 ^ [0, oo) such that e Cl ; 
(II) (j)'^ , V'^'' : M ^ M such that (f)'^ , g Cl ; 

(III) lim,^+ooC = lim,^_oo^/^?i"-0; 

(IV) The hmits 

hm (t>'^, hm Vn. 1™ Hm 



+ — ool-|-a;' x-*+ca\ — x 

exist and are bounded. 

Remark 9 Once the limits in (III) are required to exist, their vanishing is 
necessary for the initial data to satisfy M™ < oo. For the same reason we must 
have lima;^±oo J/(0,x) = and so, by the second equation in H3.6|l . 

Un{0,x)^ f [{pn+.fr,x)iO,x')-nix')]dx', 



whence (|3.9|) is necessary for having lim^^^^+oo C^(0,a;) = 0. Note also that 
Ur]{0,x) is compactly supported. 
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Besides initial conditions, uniqueness of solutions also requires data at past 
null infinity, which we take as 

(V) (/-n , 1/;^ : [0, cx)) ^ M such that 0n , V^rV G Cl n . 

Definition 2 A data set ji/jp^^n :'0n } said to he admissible if it sat- 

isfies (I)-(V) and the conditions 

(VI) hm C(x) = 0-(O), hm V^"(x)=V^(0), 

(VII) hm ^ = -0-'(O), hm ^^=^-'(0). 
Our main result is the following theorem. 

Theorem 2 Let an admissible data set {/n",<?!'n ,^/'^,'^n ,V''n } ^£ given such that 
Mq"<oo. Then there exists a unique (/n,[^nj</'n:'0n) gCoo global solution of 
such that 

(t) fn{0,x,p)^f^-{x,p), (0n,V'n)(O,x) = (C,V^'^)(x), VxeM,peM2; 

fiij lima;^_oo0n(T,a;) = 0n('r), hmx^+oo-iAnCTia;) = -011 (t), Vr^O. 

The solution of Theorem 1 is not always isolated. This depends on the choice 
of the data, according to the following theorem. 

Theorem 3 For the solution of Theorem\^we have 

Mo(r) = Af^"+ r [(0-)^ + (V^+)'l (r')dr'. (3.12) 

Moreover the solution satisfies and thus it is isolated, if and only if 

</,-(t)=V+(t) = 0,, for allT^O. 

We split the proof of the above results in five steps. In Section ITTI we de- 
rive an integral representation formula for the fields by integrating along the 
characteristics of the operators dri{lix)~^^dx in the (r,a;)-plane. The proof 
of Theorem |31 relies only upon this formula and is presented in Section 13.21 In 
Section ITSl we estimate local solutions in L°°. Unlike the case of the Cauchy 
problem considered in |13|, in the present situation we have to consider sepa- 
rately four different regions of the plane. In each of these regions, a uniform 
bound on the fields is obtained by integrating the energy identity over suitable 
subsets of the (r,a;)-plane. In Section [3.41 we formulate a local existence and 
uniqueness theorem together with a continuation criterium which states the a 
local solution can blow-up in finite time only if the momentum support of the 
particle density becomes unbounded. The proof of this result can be obtained 
by standard methods and is therefore omitted. In Section 13.51 the momentum 
support of the particle density is estimated, thereby completing the proof of 
Theorem H 
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3.1 A representation formula for solutions in the class Co 



Recall that Pq = + |pp +Pix. Since Po>0 in the support of /n under the 
assumption (|3.11|) . the Vlasov equation for solutions in the class Coo can be 
rewritten as 

Pi. . , 1 



drfn + —d^fn + — W^ + \p\''Un + (0n - i'n)p2]dp, f 
Po Po 



+ -[MV^TW~Pi) + MV^TW+Pi)]dpJ = 0. (3.13) 
Po 

The solution of (|3.13|) with initial datum f™ is given by fr]{T,x,p) = 
/p"(X(0),P(0)), where {X{s),P{s)) is the solution of the characteristic system 

i=-, Pi = - [v/1 + \p\^Un + (0n - V'n)p2], (3.14) 
Po Po 

I 



P2^-[MV^TW~Pi) + MV^Mp\^+Pi)] (3.15) 
Po 

subject to the condition {X{t),P{t)) ~ {x,p). It follows by Lemma [3 that the 
subset of the (r,x)-plane defined by 

is a vacuum region, i.e., /n — 0, for {t,x,p) gV xM?. 

Next we derive an integral representation formula for solutions of (|3.7() hav- 
ing regularity, initial data and decay as stated in Theorem |21 To this purpose 
it is convenient to introduce the auxiliary functions 

9+{z)^—;: — , g-{z) = -g+{z), z>0. 
2z 

Note that 

, , , l + z2 ±1 

Whence g± are invertible functions and we have 

g^^{z) = ^/l + z^±z, zeR. 
Consider the following subsets of the half-plane t ^ 0: 

C!i = {{t,x) : < t < mm{gZ\x),g^\x))}, ^2 = {{t,x) : gZ\x) ^ r < g+\x)}, 
fls = {(r,x) : T ^ ma.x{gZ^ix),g^^{x))}, fl^ = {(T,a;) : g+'^{x) ^ t < gZ^{x)}. 



Proposition 5 Let j'^eC^ be given such that — for {t,x)(zV and let 
^™,'!/'n')'/'n )V'n satisfying (II)-(VII) he given. There exists a unique solution 
((/)n,V'n) of ^Ijjj satisfying 
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(A) (.^n,^^n)|t=o = (C,V'^"), 

(B) lim^^_<x>?!'n(T,a;) = (j!)n(T), lim^^+oo lAnlTja;) = V'n (r), /or allr'^Q. 
Moreover the following holds: 

(i) ^n{T,x) is given by 

^n(r,a;) = (Co5+)(5;i(x)-r) 

— T M9+\y)-9-+\x)+T,y)dy, (3.16) 

^Jg+{a-^(x)-T) 

for (T,a;) eOiU02 and by 

(l>n{T,x) = 0n {r-g+^{x)) 

~l I^Jni9+\y)-9+\x)+T,y)dy, (3.17) 

for {t,x) gSIs[JQ,4. 
(a) il)n{T,x) is given by 

Mr,x) = {i^^og_){gZ\x)-T) 

I fg-{gZ^{x)-'r) 
-^J jn{9Z\y)-9Z\x) + T,y)dy, (3.18) 

for {t,x) gSIi[JQ,4 and by 

^n{T,x) = ip^ (r-gZ^ix)) 
1 f~^°° 

-^J Jn{9Z\y)-gZ\x)+T,y)dy, (3.19) 

for O2UO3. 
(Hi) In the vacuum region V we have 



(j)n{T,x)-- 



g+) (5; 1 (x) - t) , for {t,x) e (Oi U O2) n V, 
(T-5+i(a;)) /or (T,a;) e (03U04)nV, 



,1, (^^^-j i^n°9-){9Z\x)-T), /or (r, a;) e(0iUQ4)nV, 
Vn(T,x) /or(r,ar)e(02Ufi3)nV, 



fiuj The following limits 

lim 0n(''";a;) = 0, lim tpn{T,x) = 0, 

X — >-\-oc X — * — 00 

lim 0n(T,a;) = (/)n(T), lim Vn(T,a;) = ^^(r), 

x— ^ — 00 a:— ^+00 

are attained uniformly in compact subsets ofr^O. 
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Proof: We give the details of the proof for only, the argument for ipr] being 
the same. Let us show first that the function (pn defined by H3. 1 6(1 - (|3. 17(1 is . 
Owing to jp(r,a;) = 0, for Vl + a;^ '^Rq + ^t, it is easy to see that the integrals 
in the right hand sides of H3.16|l - H3.17|) are extended over compact sets and so 
(pf] is in the interior of each of the two regions ^liU^l2 and 173Uri4. The 
border between these regions is the curve T = g^^{x) and we have 

1 /"^ 

lim 0n(T,a;)= lim (j)'^{x)-- / jni9^^{y),y)dy, 



X — > — oo 



If 

lim (/)n(T,a::) = 0n(O)-- / jn{9^^{y),y)dy- 

Hence the continuity of (f>n follows by (VI) . We emphasize that the integral in 
the right hand side of the last two equations is finite, as jn{9+^{y)iy) = 0j for 
y^g-{2Ro). Taking a r-derivative in l|3.16|l and (|3.17l) we obtain 

drc^n = (r^) i9+\^)-r)-l9'A9A-)-r)fn{0,9+{9^\^)-r)) 

dTin{9+^{y)~9+^i.x)+T,y)dy, 

9+ {9+^(x)-r) 

for (T,a:) G ill U and 

1 

9r0n ='/'n'(T--5+^(a;)) - 2 j ^rjn{9+^{y)~ 9+^{x)+T,y) dy, 

for (r, a;) S ria U . Using again the support property of and (VII), we find 
that lim^ t^r^n = l™.^ c^r^n- In the same way one can show 

T >^x ) T \X ) 

that dx(t)r\ is continuous. A straightforward computation reveals that (pc] solves 
the first equation in (|3.7|l with the initial-decay conditions 0n (0,a;) = 0^ (cc), 
\\mx^-ao4'r\(.T,x) = 4>'^{T). Let us show now that any solution to the latter 
problem has the form H3.16|l - H3.17|l . thereby proving the uniqueness part of the 
theorem. Let (/)n be a solution of the first equation in H3.7|l with the initial- 
decay conditions (pni^T^) = (j)™{x), \\mx^-ao4'r\{T,x) = 4>'i^(T). It follows that 
the function h{T,z) = (/>n(T,.g+(z)), z>0, satisfies the problem 

{dr+d,)h^-^fi, /i(0,z) = (Cog+)(z), hm /i(t,z) = 0-(t), (3.20) 

Z 2^0 + 

where h{t,z) — j^{T,g^{z)). Integrating along characteristics one finds that the 
solution of H3.20|l is 

h{T,z) = {(f>'^og+){z-T)-- I {fi{s,-)g'+){z + s-T)ds, 







for T < z and 



h{T,z) = (f)f^{T-z)-^ I {fi{s,-)g'_^){z + s-T)ds, 
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for r ^ z. As h{t,z) = j'^{T,g+{z)) and by a change of variable in the integrals 
we obtain 

1 r9+(^) 

h{T,z)^{4)'^og+){z-T)-- / j^{g~'^{y)+T~z,y)dy, 
for T < z and 

I rg+{z) 

h{T,z)^cj)-{T-z)-- J j^{g~\y) + T-z,y)dy, 

for t'^z. Since 0n(''',a^) = /i(T,<?+^(a:)), this completes the proof of the repre- 
sentation formula for (/>n. Next we prove (iii). Since jp(r,a;) = 0, for vT+a?^ 
-Ro + we have 

Jn(3;'(y)-.9;'(-^)+r,y)=0 (3.21) 
for \gZ^i^y)l^RQ + \T — \g'^{x). Hence for y^x, (|3.21|l is satisfied if 

i.e., if Vl + a;2 = i(5Zi(a;)+5;^(x))^i?o + 5r. Using this in l|TTB|) - (|XT7|) con- 
cludes the proof of the claim for 0n- By (iii) and the condition (III) on the 
initial data, the proof of (iv) is straightforward. □ 

Remark 10 When 0^ =1/;^ =0, i.e., when — according to Theorem |31 — the so- 
lution is isolated, the formulas H3.17|l . H3.19|l imply that the field in the region 
r^a is given by the retarded solution of the equations. 

3.2 Proof of Theorem H 

In this subsection we prove Theorem 13 By (iv) of Proposition [3 the left hand 
side of (|3.8(l equals 



">-)%(V.+)' 



which vanishes for all Ti,r2^0 if and only if 0^ =''/'n =0- The analogue of 
(|2.18|l in one spatial dimension is 



mo(r,r)=TOo(0,r)-t- / [pn(T',-r)-pn(r',r)] dr'. 
Jo 

For r sufficientely large, the integral in the right hand side equals 



Letting r-^-^-oo and using again (iv) of Proposition |5l concludes the proof of 
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3.3 Bounds on the fields 



In this subsection we prove that the fields are uniformly bounded. Let 
(/n,J7nj^n)'0n) be a local solution of (|3.5I) " (|3.7I) which match the initial- 
decay data. Then 

{pn+jn^){T,x')dx'- n{x)dx 

-OO J —OQ 

^ No{T) + \\n\\L^=N}^ + C^C Vre[0,T), 

whence ||C/n(''')||oo ^C, for all re [0,r). Here and in the following we denote 
by C a generic positive constant which may change from line to line. The 
estimation of the fields 0n:'0n is not so simple. We shall need the following 

Lemma 3 

(en ± pn ) (r, x) (r, x), V (r, x) G [0, oo) x M. 

Proof: From one hand 



\P2\ 



-rfndp; 



on the other hand 

en 



±Pn= J (v/T+H^±pi)/ndp+^C/^ + </>^+Vn±'/'nT^^ 



The lower bound 



/y— y2±„ 1 + \p\^-pI ^ 1+pI ^ \P2\ 



yields the result. 



□ 



Proposition 6 Let (/n,C/n,0n:V'n) GCt be a local solution of f.y..5|) - f^r^ with 
admissible initial-decay data such that M™<oo. Then 

l!'/'n(T)||oo + ||V'n(T)||oo5;C, VTe[0,r). 

Proof: By Lemma O and (|3.16|) - (|3.19|l we have 

|(/>n(r,x)|<C+i / [en ~Pn) {9+' iy)-9+\x)+T,y)dy, (3.22) 

for (TjX) e J7i U il2 and 

\^n{r,x)\^C+^J {en-pn){g+\y)-g+\x)+T,y)dy, (3.23) 
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for (tjO::) £ uri4. At the same fashion, 



1 



\Mr,^)\^C+-J ien+Pn){gZ\y)~gZH^)+r,y)dy, (3.24) 
for {t,x) GfliUil4 and 

1 

\Mr,x)\^C+-J {en+Pn){9Z\y)-gZ\x)+T,y)dy, (3.25) 

for {t,x) S U fis. Thus it is enough to estimate the integrals in the right hand 
sides of (|3.22|I - H3.25|I . For {tq,xo) G fii, we integrate the energy identity H3.10|l 
in the interior of the past hght cone with vertex on (tq,xq), namely the region 
of the (T,a;)-plane so defined: 

A (to , xo ) = { (r, x) : ^ T ^ min (.g^ ^ (x) + tq - g+ ^a-Q ) , gZ ^ (x) + tq - .gZ ^ (a;o )) } • 
So doing we obtain 

{en-pn){9+^{x)-9+^{xo) + To,x) dx 

a- (9l^(^o)-ro) 

ien+Pn){gZ^ix)-gZ^{xo) + To,x) dx 

Xo 

a- {aZ''{xo)-'ro) 

{en+pnx){0,x)dxs^MiP, V(to,xo) € fli. (3.26) 

a+{a+^{xo)~ro) 

For {tq,xq) e ri2, we integrate (|3.1U|I in the interior of the truncated cone 
A(i?)(ro,a;o) = A{To,XQ)n {{t,x) : g+ (g:^^ (xq) - tq) ^x ^ R} , R>xo. 
So doing we get 

/ _ {en-pn){9+^{x)-g+^{xa)+To,x) dx 

a+{a+^ (.xo)-To) 

+ {en+pn)i9-\x)-9_\xo)+To,x)dx (3.27) 

Jxo 

rR raZ\R)+ro-aZ\xo) 

= / {en+pnx){0,x)dx- / pn(r,i?)dr. 

Ja+{a+^{xa)-To) Jo 

By (iv) of Proposition |S1 we have 

raZ^ {R)+Ta-aZ^ (xo) pTo-gZ'-ixa) ^ 

— pnir, R)dT — > / i'^n) i''')d'T, asi?— >oo. 

Jo Jo 
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Whence in the limit R^oo, (|3.27|) entails 

/>oo 

+ / (en+pn) (gZ^(a;)-.gr^(xo)+To,x) 

J xa 

s=:a4" + ||V^+|U2, V(To,xo)er!2. (3.28) 

A similar argument can be applied in the regions f^a and f24. Precisely, for 
(tojXo) € r^s the energy identity (|3.1U|) is integrated in the truncated cone 

A(-i?,fl)(T'o,a;o) = A(ro,a;o)n{(T,a;) : |a;-a;o| ^ i?}, i?>0, 

whereas for (tqjXo) G we choose 

A(__R) (tq , Xo) = A (tq , Xo) n { (r, x) : -i? < a; ^ 5_ (yZ^aio) - To) } , i? > -xo . 

So doing and letting i? — > cx) we arrive to the estimates 

{en~Pn){9+^{x)~g+^{x(i) + TQ,x) dx 



OO 

+ 00 



+ / (en+Pn)(5-\a;)-5_^(xo) + To,x) dx 

Jxo 

^M^- + U-\\l2 + \\^P+\\l2, V(To,xo)er!3, (3.29) 

(en-pn) {g+^ (x) - g:^_^ {xq) + to,x) dx 
r9- (9Z^(a;o)-To) 

{en+Pn){gZ^{x)-gZ^{xo) + TQ,x) dx 

s^M^" + ||0-||i2, V(ro,xo)efi4. (3.30) 

Combining l|X ^ -l|On |l with ijX^ - ljX^ conculdes the proof of the 

proposition. □ 

3.4 A local existence theorem 

The existence of solutions in the class Ct, for some r>0, can be proved by 
standard methods. Moreover, the solution can be continued in a larger time 
interval as long as the momentum support of the particle density is bounded. 
Precisely, we have 

Theorem 4 Let an admissible data set {f^,(j)^^,tp}^,(j)^,ip^} be given. There 
exists T>0 and a unique (/n,C/n,0n,''/'n) SCt solution of j5'.5|) - P77| ) such that 
the conditions (i)-(ii) in the statement of Theorem\^are attained. Let Tmax be 
the maximal time of existence; then 

Vn (Tlnax) < 00 =4> Tmax = OO . 
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Given the representation formulas for the fields established in Section ITTl the 
proof of the above theorem proceeds exactly as in the case of the Cauchy prob- 
lem. The argument is based upon an iteration scheme first introduced in pp. 
We refer to for details, see also ^7]. As a matter of fact, the proof in one 
spatial dimension is much easier and is essentialy given in 13', although the 
result is not explicitely stated. 

3.5 Bound on the momentum support and proof of The- 
orem [21 

Let (/n,C/nj'/'n,V'n) GCt be a local solution which match the initial-decay data. 
The uniform bounds on the fields proved in Section 13.31 lead to an estimate on 
the fimction H3.11|l . To see this, let 

sup (||C/n(T)|loo + ||0n(r)||oo + ||0n(T)|loo) 

re[0,T) 

and define 

Q(r)= sup Wl + bP + ^v/T 

+ a;2, {x,p) e supp/n(s) > >Vn{T)- 
se[o,T) J 

Along solutions of (|3.14() - (|3.15() we have 

ds\ ) Vo Po Po 

\ Po po J 

Integrating the previous inequality in the interval [0,t], t<T, we obtain the 
estimate Q{t)^C{1 + t), for all tS[0,T). If the maximal time of existence 
Imax were finite, the latter estimate would imply VniTmax) ^C{l + Tjna,x) <oo, 
which is a contradiction to Theorem ^ Thus T„iax = oo and this concludes the 
proof of Theorem [5] 

Acknov^rledgments: The author acknowledges support by the FCT, Portugal 
(contract SFRH/BDP/21001/2004). 
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